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November 1984 -April 1985 
MSG 1419* 

All of the time and effort during this report period was directed toward 
acquiring, reducing and analysing hot-wire anemometer data. This experimental 
study included twelve combinations of chord Reynolds number, angle of attack, 
and freestreao disturbance environment using the NACA 663-OI8 airfoil. 

This research has as Its objective the detailed documentation of the 
structure and behavior of the separation bubble including transition and the 
redeveloping boundary layer after reattachment over an airfoil at low Reynolds 
numbers. The intent of this work is to further the understanding of the 
complex flow phenomena so that analytic methods for predicting their forma tf on 
and development can be Improved. These analytic techniques have applications 
in the design and performance prediction of airfoils operating in the lew 
Reynolds number flight regime. 

CALCULATION OF PARAMETERS 

The primary results of this Investigation are the various flow field 
parameters that were calculated from the basic data. A majority of these 
... .parameters, were, .local, variables . calculated at the points of separation, 
transition, and reattachmont. Since the locations of these points were 
altered slightly os additional data was taken, a computer code was developed 
to handle the large number of tedious calculations. Although this calculation 
.scheme was generic in nature, the vast similarities that existed in the data 
allowed for accurate and consistent computation of the parameters. Some 

precautions had to be taken however, and those will be discussed in more 
detail. 

* NASA Technical Monitor for this Grant Is 


Mr. Den H. Scmsrs, l.'ASA Langley Research Center, 
Hampton, Virginia 23655. 
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Host of the parameters Included In Tables I and II arc derived from the 
basic boundary layer variables and require little additional discussion once 
the definition Is provided. The definitions of the parameters can be found In 
the list of nomenclature or on Figure 1. In the case of the two angular 

parameters, however, further discussion of the definition and calculation 
technique Is required. 

Separation Bubble locations 

Before any parameters could bo calculated, the locations of flow 
separation, transition, and reattachment had to be determined. These 
important locations, which collectively describe the basic character of all 
laminar separation bubbles, ware determined through a careful examination of 
the static pressure and hot-wire anemomatry data. Since static pressure data 
was cvallable at closer Intervals al eng the airfoil uppar surface. It was used 
as the primary indicator of bubble locations. The hot-wire profiles ware used 
to confirm the pressure data, and In the case of separation, ware sometimes 
used to further pinpoint the actual location. In general, the pressure and 
hot-vrlre data shewed good agreement, and together provided sufficient 
fnfc.raaticn for. determining detailed bubble locations. 

In order to determine separation, transition, and reattachsont locations 
in a consistent manner. It was sometimes necessary to choose locations between 
adjacent pressure taps. In the case of separation, this type of location 
would later permit accurate calculation of the local pressure and velocity 
gradients. The locations of transition and reattechment were taken as those 
points at which the pressure distribution exhibited sharp discontinuities. 
These discontinuities were sometimes masked because of the finite distances 
between taps. Thus, interpolation between points was necessary to simulate 
the linear pressure increase that usually exists between transition and 
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reattachment. In general, this process allowed for greater consistency In 
locating the various flow phenomena associated with the separation bubble. 

Further justification for assigning locations between taps came from the 
need for defining uncertainty bands. In most cases. It was obvious that the 
Important bubble locations occurred between a pair of adjacent pressure taps. 
When this occurred, central locations between the taps were chosen so that 
uncertainty bands extending across the two taps could be established. 
Although the exact locations of these various flow phenomena can be disputed 
to some degree, the locations determined represent mean positions within the 
bands of uncertainty. These uncertainties only apply to the pressure 
distributions obtained in this Investigation which were shown to be very 
repeatable. Long term repeatlblllty tests, however, might reveal slight 
deviations which would further extend the uncertainty bands. In general, 
small variations in separation, transition, and reattachment locations did not 
significantly affect the final results. The bubble locations associated' with 
each of the experimental cases are summarized in Table III. 

Application of Definitions 

• Although, -.the., definitions of the various flowficld parameters are 
relatively straightforward, their application to a discrete set of data is not 
necessarily trivial. On an airfoil, for example, all length dimensions should 
be maasured along the curved surface. The origin of this coordinate is 
usually taken as the point of flow stagnation. In this Investigation, the arc 
lengths along the airfoil surface were approximated by summing the. line 
segment contributions calculated from the coordinate points. Fortunately, 
additional coordinate points were available near the loading edge so that the 
linear approximation was very good. . . 

Another potential problem arose in choosing suitable values of the 









external velocities. Because of the effects of hot-wire probe orientation, 
the pressure data was used to obtain all external velocity magnitudes. Since 
the locations of separation and transition, as determined by the process 
described above, sometimes foil on gradient regions of the pressure 
distributions, the associated external velocities were measured at the extreme 
ends of the pressure plateaus Instead of at the actual locations. Although a 
relatively minor point, this distinction was necessary to combine the effects 
of consistent bubble locations with accurate parameter calculations. 

Because It was necessary to survey the entire boundary layer over the 
airfoil. It was not possible to take data at close enough Intervals to always 
correspond to the exact locations of separation, transition, and reattachment. 
For this reason, an interpolation scheme was employed to calculate values at 
the paints of interest. To determine the integrated thicknesses at 
separation, for example. It was sometimes necessary to interpolate between an 
upstream and a downstream value. Fortunately, In most cases, data was 
available at stations at or very near the desired location so -that the 
resultant Interpolated values were very reliable. 

• definition of Angular Parameters 

The turbulent spreading angle ax Is a parameter suggested by the 
simplified bubble model shown In Figure 1. In attempting to approximate the 
magnitude of this angle, the geometric simplification shown In Figure 2a was 
employed. Hare, tha "origin" of the turbulence was taken to be the point of 
maximum turbulence Intensity at transition. This point generally fell at or 
very near tha vertical center of the free shear layer. From this point, the 
turbulence via s assumed to spread between two linear boundaries extending 
outward to the boundary layer edge and the airfoil surface at reattachment. 
With this simplification, calculation of became trivial once the •'points of 
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transition and resttachment were determined. 

The final parameter which . requires additional discussion Is the angle y 
formed between the airfoil surface and the separation streamline. Because of 
the problems associated with the flow visualization data, a direct measurement 
of this angle was basically Impossible. As such. It was necessary to use 

hot-wire velocity profiles to estimate this Important parameter. This 
procedure Involved determining the heights of the recirculation region at 
three or more stations just downstream of separation. These heights were 
defined as those distances over which the velocity ratios U/Uf S were very 
small, and the profiles had extremely large slopes. These heights were 
plotted versus the surface arc length at which the corresponding profiles were 
taken. This particular definition of a separation angle on a curved surface 
Is similar to that proposed by Dobblnga et al [13. Fortunately, the points 
corresponding to the laminar free shear layer generally fell on a single line 
which could be extrapolated back to form an angle. The separation angle so 
defined is shown schematically In Figure 2b. 

DISCUSSION OF RESULTS • 

During the data acquisition phase of this investigation, the separation 
bubble flcwfield over the airfoil was surveyed for twelve different 
conditions. These conditions ware chosen in such a way that the effects of 
Reynolds number, angle of attack, and disturbance environment could be 
isolated. Along with indicating such behavioral trends, the data also 
provided valuable insight into the structure of the separation bubble. In 
addition, the development of the turbulent boundary layer downstream of 
reattachment vfas Investigated. 

The Important separation bubble parameters are compiled together in 
Tables I and II. Those values were calculated from the basic data once the 









positions of separation, transition, and reattachment were determined. A list 
of the actual bubble locations used In the analysis Is Included In Table III. 
As illustrated In Figure I, the length dimensions, lj, l2» and lg used 
throughout this Investigation represent distances measured along the airfoil 
surface. Definitions of the various parameters are Included in the list of 
nomenclature. The more Involved calculation schemes used to determine the 
angular parameters aj and y were described earlier. 

General Separation Bubble Characteristics 

As shown In Figures 3 and 4, the chord Reynolds number has an Important 
effect on the separation bubble flowfleld. These figures show the effects of 
Reynolds number for the one and zero flow restrictor cases respectively. As 
Figure 3 Indicates, the length of the separation bubble decreases as the chord 
Reynolds number Is Increased. As the tunnel speed Is elevated to achieve 
higher Reynolds number testing conditions, the rate at which small 
disturbances are being amplified In the unstable laminar shear layer is 
Increased. This causes forward movement In the point of transition' which 
subsequently reduces the overall length of the separation bubble. 

A$ the bubble . dltninlshes, In length, the pressure distribution becomes 
less "distorted", and a higher suction peak is attained. This implies- that 
changes In the separation bubble significantly affect the entire leading edge 
flow. As a result, the upstream pressure distribution is highly dependent on 
the chord Reynolds number. As shown In Figure 4, however, the effect of the 
bubble on the downstream flowfleld appears to be minimal. 

The effect of chord Reynolds number on total bubble length Is summarized 
In Figure S. In this figure, the bubble lengths lg listed In Table I are 
plotted versus Reynolds number for all of the 12 degree angle of attack cases. 
As Indicated, the bubble more than doubles In length as the chord' Reynolds 
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number is decreased from R c ■ 140,000 to R c • 50 , 000 . Each of the different 
flow restrictor cases exhibit the sarse trend with Reynolds number although the 
magnitudes of lg are shifted. This can be seen by the dotted line which 
connects the three cases at R c ■ 140 , 000 . The shift magnitude of lg 
indicates that bubble length is also a direct function of the disturbance 
environment. Although the level of freestream turbulence Increased as the 
chord Reynolds number (tunnel speed) was Increased, the latter effect appears 
to predominate. It is, however, impossible to totally uncouple the two 
effects when changes in tunnel speed are used to alter the chord Reynolds 
number. 

Although changes in Reynolds number could not be achieved without slight 
changes in the disturbance environment for the given model, it was possible to 
isolate the effects of freestream turbulence. This was made possible through 
the introduction of flow restrictors which have been shown to alter the 
testing environment [2]. Pressure distributions obtained at various levels of 
freestream turbulence are shown in Figures 6 and 7. As the turbulence level 
is Increased, the bubble Is reduced In length, and the suction peak grows in 
absolute magnitude. This well known phenomenon so closely resembles those 
effect obtained by increasing the chord Reynolds number that the two- are often 
equated. This has given rise to the use of "effective" Reynolds numbers 
when dealing with various disturbance environments or when discussing the 
effects of different types of surface roughness [3]. 

The correlation between Increases in freestream turbulence and Increases 
In Reynolds number is shown dramatically in Figure 8. Except for the slight 
deviation In the suction pecks, these distributions arc basically identical. 
Thus, the case of R c ° 140,000, a 3 12°, 2 flow restrictors can be said to 
have an "effective" chord Reynolds number of 200,000 with 0 flow restrictors. 
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Starting frost ths base condition (R c » 140,000, OFR), an increase in 
turbulence intensity of 6122 or an Increase in Reynolds number of 43? is 
required to produce the common pressure distribution. Associated with the 
Reynolds number Increase was a 42? Increase in the freestream turbulence 
Intensity. Given these Incremental magnitudes. It appears that. In general, 
the small rise In turbulence level associated with increases In tunnel speed 
contributes relatively little to the overall effect of increasing the chord 
Reynolds number. Thus, for a given tunnel configuration, changes In speed 
predominate over accompanying changes in turbulence levels as the driving 
mechanism which affects the separation bubble. 

Along with chord Reynolds number and freestream turbulence intensity, the 
separation buoble flowfleld is also affected by changes In angle of attack. 
These effects are shown In Figures 9-11 for three different testing 
conditions. As the angle of attack Is increased from 8 to 10 degrees, the 
point of laminar • separation moves forward from approximately 3.7? to 
approximately 2.8? X/C, but there is no significant change in the length of 
the bubble. At 12 degrees angle of attack, the bubble has again moved forward 
but has now Increased in length. As a given angle of attack, the separation 
point reami ns essentially unchanged over the entire range of Reynolds number. 
The effect of angle of attack on the bubble length lb is summarized in Figure 
12 . 

In addition to the overall length, the thickness of the separation bubble 
Is a parameter that Is significantly affected by the various testing 
conditions. For the sake of- comparision, the height of the recirculation 
region at transition (Hg)x was determined from each set of velocity profiles. 
These thicknesses are plotted versus chord Reynolds number In Figure 13. This 
figure Indicates that a sharp reduction in bubble thickness occurs '-as the 



C O 


e. 


o . Q 






10 


chord Reynolds number Is Increased. The bubble . at R c =» 50,000 Is almost 3 

times thicker at transition than the corresponding bubble at R c ■ 140,000. 
This seems to Imply that the bubble thickness is closely related to the total 
bubble length. This almost linear relationship is plotted In Figure 14. The 
limited number of points available in Table I also Indicate that the bubble 
thickness Increases as the angle of attack Is Increased. This phenomenon was 
confirmed by flow visualization data. 

In addition to determining the overall separation bubble characteristics, 
the recent work has focused on documenting the various structural components 
comprising the bubble. This documentation Involved a detailed investigation 
of the various flow phenomena associated with the bubble, as well as, a study 
of the rede\ loping turbulent boundary layer downstream of reattachment. Once 
this documentation was completed, the largo amount of experimental data was 
analyzed In terms of existing physical and mathematical models of the flew 
field. The results of this analysis Indicated that further work Is needed to 
adequately model the flow over an airfoil at low Reynolds numbers. It 
appears, however, that this work, which Is described In detail in Reference 
*• [2], .has'- : provided -the groundwork for future Improvements In the design and 
performance prediction of low Reynolds number airfoils. 
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NOMENCLATURE 


C Airfoil chord, mm 

C p Pressure coefficient, C p » (?^ - Pf s )/Qf s 
PR Flow Restrictor 
Hy Transition height, m 

H12 Shape factor eaual to the boundary layer displacement thickness 
divided by the momentum thickness 

H32 Shape factor equal to the boundary layer eneroy thickness divided by 
the momentum thickness 
ig Total bubble length, r.tn 

*1 Bubble length from separation to transition 

z? 3 ubble length from transition to reattachment 


R $2 



Reynolds number based on momentum thickness, U / P 
Reynolds number based on laminar length, T J 1^ /5 / P 

Reattachment location 

Reynolds number based on displacement thickness, U b^P/P 
Chord Reynolds number 

Reynolds number based on surface arc length 
Reynolds number based on total bubble length 


'S' laminar separation Tocation - ‘ Surface arc length coordinate 


S’ Turbulent separation location 

T Location of approximate end of transition 

U Velocity 


U' Fluctuating velocity component 

; Jf S Freestream velocity 

X/C Nondimansional distance along chord 



Greek Symbols 
a Angle of attack 
y Separation angle 

Turbulent spreading angle 
^ Boundary fayer displacement thickness 
$3 Boundary layer energy thickness 
S 2 Boundary layer momentum thickness 
v Kinematic viscosity 
P Density 

Subscripts 
fs Freestream 
R Re attachment 
S Separation 
T Transition 
ext External 
atm Atmosphere 
v . Uncertainty magnitude 
, E Bubble, 

S' Turbulent Separation 







TABLE I 

IMPORTANT FLOW FIELD PARAMETERS 
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R c *ICOX 
0-12' 
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Rr»20OX 
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OFR 
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o*12*- 
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1 

t s (S/C Si 

13.6 

12.3 

9.0 

22.0 : 

12.7 

il {S/C S} 

7.1 

6.4 

«.3 

12.9 

8.3 

12 (S/C I) 

6.5 
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4.2 

9.1 \ 

4.4 

(h 3 ) t <:n 
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0.29 
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0.46 
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17.0 
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17.1 
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25.4 
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9.2 

Up, (rt/s) 
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20.3 

4.5 

7.4 

(ii)z (SC) 

C.lo 

0.13 

0.14 

0.25 

C.23 

f5 2 } a CCl 

0.53 

0.66 

0.35 

1.37 

0.45 

(527s i *C) 
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0.059 
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0.20 
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0.67 

11.7 
IX. 7 
9.5 
0.19 
0.41 
0.052 
0.23 
3.50 
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a c «ioox 

o-12* 

1FR 


14.2 
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5.4 
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12.0 

12.1 

9.5 
0.23 
0.75 
0.043 
0.39 
4.91 
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o”12* 

1FR 


10.3 
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0.37 
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17.3 
17.2 
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0 . 1 <» 

0.«4 

0.044 

0.26 

4.28 

1.71 


R c -5« 

o*10* 

2FSJ 


15.5 

10.0 

5.5 
0.67 
C-.86 
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0.26 
0.31 
0.079 
0.35 
3.26 
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1.29 
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1.34 
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2.20 
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1.00 
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4.21 

5.06 

1.06 

2.23 
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509.0 

229.0 

346.0 

109.0 

105.0 

150.0 

55.0 

07.0 

29.0 

25.5 

27.5 

35.5 

30.0 

e.o 

6.5 

6.0 

9.0 

4.0 

2.56 

2.79 

2. 69 

1.87 

2.40 

1.53 

1.52 

1.53 

1.63 

1.55 


1.64 

1.61 

2.42 

1.34 

1.02 

2.36 

343.0 

72.0 

26.5 

6.5 

2.10 
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r c *1ooj: 

n* 10* 

IFK 

R c «loor. 

o*12‘ 

1FR 

R c *140r 

o t 12* 

IFR 

R c -5rac 

a* 10* 

?FR 

R c * 5(X 
a*l 2* 
2FR 

1.51 

1.40 

1.50 

1.50 

1.53 

1.66 

1.61 

1.66 

1.55 

1.57 

1.06 

2.64 

2.00 

1.40 

1.73 

1.26 

1.44 

1.40 

0.003 

1.01 

1.24 

1.31 

1.07 

0.671 

n.49o 

2.50 

2.75 

3.46 

1.57 

1.71 

346.0 

424.0 

514.0 

233.0 

226.0 

95.0 

79.0 

119.0 

71.0 

67.0 

51.0 

20.0 

25.5 

30.5 
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5.0 

6.5 

6.0 

3.0 

6.6 

2.52 
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SEPARATION BUBBLE LOCATIONS 


Condition 

Separation 
(X/C *) 

Transition 
(X/C %) 

Reattachment 
(X/C %) 

R c »140,000 
Aipha-12 Ocg 
0 FR 

2.0 +/- 0.4 

8.7 +/- 0.6 

15.0 +/- 1.0 

R c «160,000 
Alpha-12 Deg 
0 FR 

2.0 +/- 0.4 

8.0 +/- 0.7 

13.7 +/- 1.0 

R--200.000 
Alpha-12 Deg 
0 FR 

2.2 +/- 0.2 

6.7 +/- 0.6 

10.7 +/- 0.8 

R c -50,000 
Alpha-12 Oeg 
1 FR 

2.0 +/- 0.4 

14.3 +/- 0.6 

23.3 +/- 1.2 

R r *8Q,00Q 
Afpha-lO Deg 
1 FR 

2.8 +/- 0.5 

10.7 +/- 0-3 

15.0+/- 0.7 

R_=80 ,000 
Alpha-12 Deg 
1 FR 

2.0 +/- 0.4 

10.7 +/- 0.8 

18.0 +•/-• 1.2 

R c «100,000 
• Atpha-1Q. Deg ■ 
1 FR ' 

2.8 .+/- 0.5 

9.3 +/- 0.6 

12.5 +/.- 0.6 

R--100 ,000 
Alpha- 12 Deg 
1 FR 

2.0 +/- 0*4 

9.3 +/- 0.6 

15.6 +/- 1.2 

R c =140,000 
Alpha-12 Oeg 
1 FR 

2.2 +/- 0.2 

7.3 +/- 0.4 

12.0 +/- 0.5 

R c -50 ,000 
Alpha-10 Deg 
2 FR 

3.4 +/- 0.4 

13.0 +/- 0.5 

18.4 +/- 1.2 

R c «50,000 
Alpha- 12 Deg 
2 FR 

2.5 +/- 0.3 

13.0 +/- 0.5 

20.7 +/- 1.2 

R_*140,000 
Alpha- 12 Deg 

2.2 +/- 0-2 

6.7 +/- 0.6 

10.7 +/- 0.8 




SEPARATION TRANSITION REATTACHMENT 


Flguro 1, Two -Dine r.sional Separation Bubble Model 
Sucsested by Sppler 
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Figure 2. Calculation Techniques for Angular Thrsunotars 
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1 Flow Restrictor 
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Figure 6 . Effect of Freestreaa Turbulence Intensity on Pressure 
Distribution, Rc = 1^0,000, a = 12 
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Figure 14. Bubble Height at Transition versus Total ' 
Bubble Length 









End of Document 


r 



